Redistributing Chern numbers of Landau subbands: tight-binding electrons under 
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We investigate the magneto-transport properties of electrons on a square lattice under a magnetic 
field with the alternate flux strength <j> ± A<f> in neighboring plaquettes. A new peculiar behavior of 
the Hall conductance has been found and is robust against weak disorder: if tf> = ^ x 2n (p and 2N 
are coprime integers) is fixed, the Chern numbers of Landau subbands will be redistributed between 
neighboring pairs and hence the total quantized Hall conductance exhibits a direct transition by 
±Ne 2 jh at critical fillings when Acj> is increased from up to a critical value A(f) c . This effect can 
be an experimental probe of the staggered-flux phase. 
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PACS numbers: 73.43.Cd, 75.47.-m, 73.20.Fz, 71.10.Hf 

Introduction. — The integer quantum Hall effect 
(IQHE) is observed in two-dimensional (2D) electron 
systems, in the presence of a strong perpendicular 
uniform magnetic field [lj. The phenomenology of 
IQHE is well explained by Landau quantization and 
disorder-induced localization. The IQHE electronic 
states exhibit interesting topological properties [1, 0] . In 
particular, each state can be labelled by a topologically 
invariant integer called the Chern number, which is the 
Hall conductance of this state, in units of e 2 /h. A state 
with a nonzero Chern number carries a Hall current and 
is necessarily extended. In solving the problem regarding 
how an IQHE state transforms into the insulating state, 
many authors carried out numerical studies of the tight- 
binding model in the presence of a uniform magnetic 
field and a random potential 0, H, @, 0] ■ An interesting 
behavior has been reported by Sheng and Weng that 
disorder will induce direct higher IQHE plateau to 
insulator transition [g|. These theoretical results are 
in good agreement with several magneto-transport 
experiments performed in Si metal-oxide-semiconductor 
field-effect transistors Q and a 2D hole system confined 
in a Ge/SiGe quantum well 

Recently, an intriguing orbital-current-carrying state 
in the square lattice, the staggered-flux phase (SFP) 
which is also known as orbital antiferromagnet or d- 
density wave (DDW) state 0, regains new atten- 
tions since unusual experimental findings in two sys- 
tems: a pseudogap phenomena (See the review by 
Timusk et al. [Ill]) in the underdoped region of high-T c 
cuprates (Io| : a hidden order in the heavy- fermion com- 
pound URu2Si2 [Hj]. However, a direct experimental ob- 
servation of the SFP order is difficult. On the other hand, 
adopting a superconducting wire network decorated with 
an array of ferromagnets, an artificial system of tight- 
binding electrons under a staggered-modulated magnetic 
field has been realized [l3], and the Little-Parks oscil- 
lations exhibit the edge part of the corresponding Hofs- 
tadter spectrum. 



Motivated by the above experimental and theoretical 
studies of the IQHE and the SFP, we investigate the 
magneto-transport properties of tight-binding electrons 
on a 2D square lattice in the presence of a staggered- 
modulated magnetic field which has the alternate flux 
strength cj> ± A<fi in neighboring plaquettes. Hofstadter's 
butterflies of similar systems have been studied numer- 
ically 
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Here, we find a new peculiar Acj)- 
dependent redistribution behavior of Chern numbers of 
Landau subbands and propose that this effect be an ex- 
perimental probe to reveal the possible internal orbital- 
antiferromagnetic structure of the SFP order in some sys- 
tems 
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1: (color online). 
4> + A(j> and <f>2 = < 



Illustration of the selected gauge. 
-A<f> 



Formulation. — The tight-binding model Hamiltonian 
is given as follows: 



H 



(i) 



where the hopping integral t is taken as the unit of en- 
ergy, Ci (cj) is an electron annihilation (creation) opera- 
tor on site i, and (ij) refers to two nearest neighboring 
sites. The magnetic flux per plaquette (the summation 
of aij along four links around a plaquette) is given by 
cf>i = 4> + A<j) and <j>2 — <t> — A<f) alternatively in neigh- 
boring plaquettes. <j> — 2W x ^ 7r (P an< ^ are coprime 
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integers) measures the uniform-flux strength, A(j> mea- 
sures the staggered- flux strength, and both are in units 
of <fio/2ir (4> = hc/e is the flux quantum), is a random 
potential uniformly distributed between [— W/2, W/2]. 

The gauge selection for this staggered-modulated mag- 
netic field is shown in Fig. [1] and the corresponding peri- 
odical boundary condition is adopted. In the absence of 
disorder (W — 0), the size of the chosen magnetic unit 
cell is 2N x 2 for <j> — ^ x 27r. Given N, the variation of 
A(j) does not change the periodicity of the system. In the 
presence of disorder (W ^ 0), the sample considered has 



the 27V x 2N unit cell for . 



2..V 



x 2tt and the periodical 



boundary condition is also adopted. 

After the numerical diagonalization of the Hamiltonian 
[Eq. JTJ], the total Hall conductance for the system can 
be calculated through the standard Kubo formula: 



ME) 



e m <Ee n >E 

(m\v x \n) (n\v y \m) 



\n)(n\v x 



) 

where Aq is the area of the system, E is the Fermi energy, 
e m and e„ are the corresponding eigenvalues of the eigen- 
states \m) and \n). The velocity operator is defined as 

i<Pi + r,i _ r 1 r _. , _r-*^<+t.*) with r = x or 

C n 



y. And we can rewrite oh as cr-n(E) 



where C m is the Chern number [2j, |3[ of the m-th sub- 
band. 



X-1 






(b) AiH> 
+1 


+1 


rk 




(c) Afl»=0.1jT 
+1 


+1 






(d) A*=0.15it 
+1 


+1 


Critical! I 




(e) A(jt=0.27i 

+1 


+1 





■4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 



-4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 

E E 
FIG. 2: (color online). The case with p = 1, N = 4 and 
W = 0. (a) The Hall conductance <th versus the Fermi energy 
E for various A0's. (b)-(e) The DOS of the cases in (a). The 
Chern numbers of subbands are shown. 

Absence of disorder, (a) An example with p — 1 and 
N = 4. — An overall picture for the Hall conductance cal- 
culated by Eq. © is shown in Fig. 02a) with the uniform- 
flux strength <f> = | x 2ir and various staggered-flux 
strengths A</> = 0, O.Itt, 0.14tt, 0.15tt, 0.2tt (only E < 
part is shown here). For the uniform case (A(j) = 0), three 
well-defined IQHE plateaus at <th = le 2 /h (I — 1, 2, 3) are 
clearly shown, corresponding to four Landau subbands 
[see the corresponding DOS in Fig. [2jb)] centered at the 
jumps of the Hall conductance at E < 0. With A0 in- 
creasing from to 0.27T, one sees a systematic evolution 



of these IQHE plateaus. At A0 = O.Itt, the / = 3 IQHE 
plateau narrows from both sides, and the I = 2 IQHE 
plateau narrows from the right side. At A<fi = 0.2ir, 
the I = 3 IQHE plateau disappears and is replaced by a 
I = — 1 IQHE plateau. A<fr — 0.157T is close to the critical 
point where the transition from the / = 3 IQHE plateau 
to the / = — 1 IQHE plateau takes place. All the values 
of the Hall conductances in Fig. [2a) fall back to zero 
when the Fermi energy approaches E = 0, due to the 
particle-hole symmetry. 

Such a behavior can be explained through the evolu- 
tion of the density of states (DOS) depicted in Fig.[^b)- 
(e). In Fig. 02b), (c) (A</> = 0, O.Itt), there are four 
separated subbands (each totally-filled subband con- 
tributes I to the filling factor) at E < with the cor- 
responding well-defined four sequential Chern integers 
{+l,+l,+l,-3}. In Fig. EDJd) (A<P = 0.15%), the up- 
per two subbands merge together and form a pseudogap 
which makes the corresponding quantization of Hall con- 
ductance not well defined. In Fig. OJe) (A</> = 0.2tt), 
the upper two subbands separate again and the sequen- 
tial Chern integers are redistributed as {+1, +1, — 3,+l}. 
We can conclude that at the critical filling v = |, there 
is a direct transition of the Hall conductance from the 
I = 3 IQHE plateau to the I = — 1 IQHE plateau when 
A<j) is increased from up to a critical value A<p c . 

(b) Redistribution behavior of Chern numbers. — Kip 
1 . we have numerically confirmed that the Hall conduc- 
tances at weaker flux strengths with N — 7-36 all ex- 
hibit the similar features as those in Fig. [2^a) . Thus in 
general, given <p = x 2 71 ", the N sequential Chern in- 
tegers {+1, . . . , +1, +1, —A + 1} will be redistributed as 
{+1, . . . , +1, — N + 1, +1} (one can also say that there is 
a transfer of integer N between the last two Chern num- 
bers) when A(f> is increased from up to a critical value 
Acf> c . Namely, there is a direct transition of the Hall con- 
ductance from the I = N — 1 IQHE plateau to the I = — 1 
IQHE plateau) at the critical filling v c = ^jf-- This is 
quite different from the global phase diagram proposed 
by Kivelson, Lee and Zhang [3] for IQHE which predicts 
that the only allowed transitions are the nearest-neighbor 
plateau-plateau transitions (/ — > / ± 1), and is also dif- 
ferent from the disorder-induced higher IQHE plateau to 
insulator transition (i — >• 0) @. 

In cases with p > 1, based on a Diophantine equa- 
tion , a pattern of subband structure has been found 
by Yang and Bhatt [7[ : the subbands away from the band 
center form groups with p subbands each, and the total 
Chern number of each group is +1, while the subbands 
near the band center form a special group with mod(N,p) 
subbands. 

At p = 3 and various N's, through numerical 
calculations, we find a more complicated redistribu- 
tion behavior of Chern numbers (Table Jj). For in- 
stance, we take the case N = 8, when A<j) varies 
from to 0.2tt, the eight sequential Chern numbers 
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TABLE I: Redistribution behavior of Chern numbers (of the 
E < subbands) for p — 3, 5, 7, 9 and various N's. Only the 
last two or three groups are shown for N > 13. 



p 


A r 




Chern numbers of Landau subbands 


3 


8 




0.037T 

O.br 

0.27T 


(_5 ii _5)(_5 ii _5)(_5 3) 
(-5,ll,-5)(-5,3,3)(-5,3) 
(3,3,-5)(-5,3,3)(-5,3) 

V I ' /V 7 7 / \ " / 

(3,3,-5)(-5,3,3)(3,-5) 


3 


10 




0.27T 


(7, -13, 7)(7, -13, 7)(7, -13, 7)(-3) 
(-3, -3, 7)(7, -3, -3)(-3, -3, 7)(-3) 


3 


11 




0.157T 


(-7, 15, -7H-7, 15, -7)(-7, 15, -7)(-7, 4) 

\ "77 /V ■ 7 " / V 7 7 / \ 7 / 

(-7,4,4)(4,4,-7)(-7,4,4)(4,-7) 


3 


13 




O.Itt 


... (9, -17, 9) (9, -17, 9)(9, -17, 9) (-4) 
. . . (-4, -4, 9)(9, -4, -4)(-4, -4, 9)(-4) 


3 


16 




O.Itt 


... (11, -21, 11)(11. -21, ll)(-5) 
...(ll,-5,-5)(-5,-5,ll)(-5) 


3 


17 




0.087T 


. . . (-11, 23, — 11)(— 11, 23, — 11)(— 11, 6) 
...(6,6,-ll)(-ll,6,6)(6,-ll) 


3 


32 




0.057T 


. . . (-21, 43, -21)(-21, 43, -21)(-21, 11) 
. . . (11, 11, -21)(-21, 11, 11)(11, -21) 


5 


12 




0.27T 


(5, 5, -19, 5, 5)(5, 5, -19, 5, 5)(5, -7) 
(5,5,-7,-7,5)(5,-7,-7,5,5)(-7,5) 


5 


13 




0.27T 


. . . (-5, -5, 21, -5, -5)(-5, -5, 8) 

\ 3 7 7 " / V ' 7 / 

...(-5,-5,8,8,-5)(-5,8,-5) 


5 


14 




0.27T 


. . . (-11, 17, -11, 17, -11K-11, 17, -11, 3) 

\ 7 3 7 ? / \ 7 5 " / 

...(-11,3,3,3,3)(3,3,-11,3) 


5 


17 




0.087T 


... (7, 7, -27, 7, 7)(7, 7, -27, 7, 7)(7, -10) 
... (7, 7, -10, -10, 7)(7, -10, -10, 7, 7)(-10, 7) 


5 


33 




0.047T 


. . . (-13, -13, 53, -13, -13)(-13, -13, 20) 
. . . (-13, -13, +20, +20, -13)(-13, 20, -13) 


7 


16 




0.087T 


. . . (-9, 23, -9, -9, -9, 23, -9)(-9, 7) 
...(-9,7, 7,-9,-9,7,7)(7,-9) 


7 


23 




0.067T 


. . . (-13, 33, -13, -13, -13, 33, -13)(-13, 10) 
. . . (-13, 10, 10, -13, -13, 10, 10) (10, -13) 


7 


30 




0.047T 


. . . (-17,43, -17, -17, -17, 43, — 17)(— 17, 13) 
. . . (-17, 13, 13, -17, -17, 13, 13)(13, -17) 


7 


31 




0.047T 


... (9, 9, 9, -53, 9, 9, 9)(9, 9, -22) 
...(9,9,-22,-22,9,9,9)(9,-22,9) 




on 




0.067T 


... (9, 9, -31, 9, 9, 9, -31, 9, 9)(9, -11) 
... (9, -11, -11, 9, 9, -11, -11, 9, 9)(-ll, 9) 


9 


22 




0.067T 


... (5, 5, 5, 5, -39, 5, 5, 5, 5)(5, 5, 5, -17) 
... (5, 5, 5, -17, -17, 5, 5, 5, 5)(5, 5, -17, 5) 


9 


29 




0.047T 


(13, 13, -45, 13, 13, 13, -45, 13, 13)(13, -16) 
(13, -16, -16, 13, 13, -16, -16, 13, 13)(-16, 13) 


9 


38 




0.037T 


(17, 17, -59, 17, 17, 17, -59, 17, 17) (17, -21) 
(17, -21, -21, 17, 17, -21, -21, 17, 17) (-21, 17) 



{(-5, +11, -5), (-5, +11, -5), (-5, +3)} (each pair of 
parentheses represents a group according to the group- 
ing scheme by Yang and Bhatt Q) are redistributed as 
{(+3, +3, -5), (-5, +3, +3), (+3, -5)}. We note that, in 
each group, one Chern number is subtracted by 8 and an- 
other neighboring Chern number is added by 8. Hence, 
the total quantized Hall conductance exhibits a direct 
transition by ±8e 2 /h at the critical fillings v c = jg, jg 
and jq. At N = 10, there is only one Chern number in 
the last group which remains unchanged when A<j) varies, 
while the other three groups change similarly as the case 
N = 8. We have numerically confirmed that for p = 5, 
7 and 9 (Table [J), there are also redistribution behaviors 
similar to p = 3. 

In general, at p > 3, there is a redistribution behav- 
ior of Chern numbers among each group [except the last 
group when mod(A^, p) — 1 and there is only one Chern 
number within it]: a transfer of integer TV between two 
neighboring Chern numbers, and the redistributed Chern 
numbers of neighboring groups away from the band cen- 
ter are antisymmetric in sequence. Hence, the total quan- 
tized Hall conductance will exhibit a direct transition by 
±Ne 2 /h at critical fillings. 

(c) Scaling properties. — Both the uniform-flux 
strength <fi = ^ x 2tt at N = 4-38 and the staggered- 
flux strength in the previous calculations are relatively 
high fields. In order to meet the requirements of possible 
experimental realization, we will reduce the flux strength 
further. The scaling behavior of A0 C with (f> is shown 
in Fig. [3l For p > 3, we concentrate on the cases where 
there are two Chern numbers in the last group and 
consider the critical filling v = -^r- For a fixed p, 
A(j) c scales with 4> almost linearly. The slopes of the 
lines with different p's differ from one another. With 
4> deceasing, all lines approach the origin point, and 
therefore A(f> c approaches zero at the weak- field limit. 
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FIG. 3: (color online). A0 C versus <f> for p = 1, 3, 5, 7, 9 and 
various iV's (at the critical filling v = J ^~ 1 ) 



Presence of disorder. — In realistic solids, impurity 
or phonon-induced disorder always has effect on the 
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FIG. 4: (color online). The case with p = 1, N = 4 and 
W — 1 (1000 random-potential configurations of the lattice 
size 8x8). <th versus E for various A0's. The inset shows 
the evolution of oh versus A(f> at v = |. 

As for p — 1 and iV = 4, we can see from Fig. |4] that 
the presence of weak disorder does not smear out the 
transition of the Hall conductance from the / = 3 IQHE 
plateau to the I = — 1 IQHE plateau near the critical 
filling v — |. 

An experimental probe of the SFP order. — As one ap- 
plication of our theoretical concern, we would take the 
staggered-flux part ±A</> as the possible internal orbital- 
antifcrromagnetic structure of the SFP order in some 
strongly-correlated electron systems [13, EH, and the 
uniform-flux part <f> as an external experimental probe. 
In Fig. at the electron filling v — ||, we fix the value 
of A<p and do such a hypothetical experiment by mea- 
suring the Hall conductance which varies with (f>. The 
two curves of the Hall conductances with A(j> = and 
A<p = 0.17T exhibit large deviations at many values of 
N/p. Therefore, such an experiment could tell us the 
characteristic differences between the states with and 
without the SFP order. 

Summary and discussion. — As for the tight-binding 
electrons on a square lattice under a staggered-modulated 
magnetic field with the strength <f> ± A0, a new peculiar 
behavior of the Hall conductances has been found: (a) 
if <fi = x 2ir (p and 2N are coprime integers) is fixed 
and A<p is increased from up to a critical value A<^ c , 
the Chern numbers of Landau subbands will be redis- 
tributed between neighboring pairs; (b) and hence the 
total quantized Hall conductance exhibits a direct tran- 
sition by ±Ne 2 /h at one or several critical fillings; (c) 
the neighboring two subbands merge together and form 
a pseudogap when A</> = Acj) c ; (d) at a fixed p, A0 C scales 
with almost linearly; (e) weak disorder does not smear 
out the peculiar transition of the Hall conductance. 

To observe this effect under an accessible uniform mag- 



netic field with N ~ 1000, we estimate that both the 
strength of disorder W and the temperature faT should 
not exceed ^ ~ ^ono ( wnere 4t is half of the original band 
width). Taking t = 0.01 eV, the temperature T should 
be lower than 0.5 K. 

As for high-T c cuprates, the strength of the internal 
staggered magnetic field is estimated to be about 0.1 
T |l7|. One may suspect that: will an internal staggered 
magnetic field survive under a stronger external uniform 
one? We believe it does, basing on the numerical re- 
sults: an external uniform magnetic field will suppress 
the cZ-wave superconducting order, while the SFP order 
(or equivalently, the DDW order) can exist or even be 
enhanced in the vortex cores [H, EH • 
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